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general  inear  complementarity  problem  has  bebn  shown  to  be  strongly  NP-Complete. 
Efficient  methods  have  been  developed  for  checking  the  convexity  of  a  piecewise 
linear  function  defined  on  a  convex  polyhedron,  in  many  parametric  cost  0-1 
integer  and  network  flow  problems,  the  number  of  sl6pe  changes  in  the  optimum 
objective  value  function  has  been  shown  to  grow  exponentially  with  the  size  of 
the  problem.  Perturbation  results  were  obtained  for  lih^ar  complementarity 
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SUJttARY  OF  RESEARCH  ACCOMPLISHMENTS. 

1.  COMPUTATIONAL  COMPLEXITY  OF  ALGORITHMS  FOR  LINEAR  COMPLEMENTARITY 
PROBLEM  (LCP).  } 

An  LCP  of  order  n  Is  a  problem  of  the  form:  find  v  e  Rn,  z  c  Rn 
satisfying 

v  -  Mz  -  q 

w  >  0,  z  >  0  (1) 

T 

wz  -  0 

where  q,  M  ere  given  matrices  of  orders  n  x  1  and  n  x  n  respectively. 

The  first  work  in  this  area  Is  a  paper  by  K.  G.  Marty  In  which  he 
constructed  a  class  of  LCPs  of  orders  n  >  2,  associated  with  lower  triangular 
P-matrices  (which  are  positive  aemideflnlte.  but  not  positive  definite), 
and  established  that  either  Lemke's  complementary  pivot  method  or  Murty's 
least  index  principal  pivoting  method  require  2n  pivot  steps  to  solve  the 
nth  problem  In  the  class,  before  termination. 

In  [2]  Y.  Fathl  established  that  these  two  algorithms  exhibit  the 
same  exponential  growth  computational  requirements  In  the  worst  case  even 
on  LCPs  associated  with  positive  definite  symmetric  matrices  M. 

In  (11]  J.  Birge  and  A.  Gana  established  that  the  Cot tle-Dan trig 
principal  pivoting  algorithm*  and  the  Van  der  Heyden  variable  dimension 
algorithm  both  exhibit  exponential  growth  computational  requlrments  In 
the  worst  case  •  on  LCPs  associated  with  P-matrlces  M. 
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2.  RESULTS  ON  LCPs  ASSOCIATED  WITH  P-MAT^C^,,  J> 

Chief,  Technical  Information  Division 

Let  H  be  a  P-matrlx  of  order  n  and  J  £{l,...n).  Consider  the 

j 

generalized  LCP  denoted  by  (q,  M,  J):  find  w  e  Rn,  z  e  Rn  satisfying 

w  -  Mz  ■  q 

T  n 
w  z  ■  0 


Wj*  *J  »  ®  for  J  i  J 
<  0  for  all  j  e  J 


(2) 


In  [1]  K.  6.  Murty  has  proved  that  (q,  H,  J)  has  a  unique  solution  for 
each  q  e  Rn,  iff  M  is  a  P-matrix. 

Assume  that  q  is  nondegenerate  in  (1)  and  that  M  is  a  P-aatrix.  In 
[1]  K.  G.  Murty  has  proved  that  given  any  +,  -  sign  vector  in  Rn,  there 
exists  a  unique  coapleaentary  basis  for  (1)  such  that  the  components  in 
the  updated  q  with  respect  to  it  correspond  to  the  specified  sign  vector. 
This  establishes  a  1-1  correspondence  between  complementary  bases  for  (1) 
and  sign  vectors  in  Rn. 

3.  COMPLEXITY  OF  THE  LCP. 

In  [5]  S.  J.  Chung  has  proved  that  the  general  LCP  is  strongly 
NP-complete. 

4.  THE  NEAREST  POINT  ALGORITHM. 


Let  A  be  a  nonsingular  square  matrix  of  order  n  and  b  a  given  point 
in  R° .  For  finding  the  nearest  point  (in  terms  of  the  usual  Euclidean 
distance)  in  the  slmplicial  cone  Pos(A)  to  b,  K.  G.  Murty  and  Y.  Fathi 
[4]  have  developed  an  efficient  algorithm.  This  algorithm  can  also  be 
used  to  solve  LCPs  associated  with  positive  definite  symmetric  matrices. 


5.  COMPUTATIONAL  COMPLEXITY  OF  PARAMETRIC  LINEAR  PROGRAMMING 


In  [3],  K.  G.  Murty  has  shown  that  the  total  number  of  slope 

• 

j 

changes  In  the  optima  objective  value  function  of  a  parametric  linear 
program  grows  exponentially  with  the  problem  size  in  the  worst  case. 

6.  EFFICIENT  BLOSSOM  ALGORITHM  FOR  EDGE  COVERING  PROBLEMS. 

Let  G  •  (N,  A)  be  a  connected  undirected  network  with  e  ■  (c^) 
as  the  vector  of  edge  costs.  A  subset  of  edges  X£^A  is  an  edge  cover 
In  G  if  for  each  i  e  II,  there  exists  at  least  one  edge  in  E  incident  “ 
at  node  i. 

In  [6] ,  K.  G.  Murty  and  C.  Perin  developed  an  efficient  blossom 
algorithm  for  the  problem  of  finding  a  cost  edge  cover  in  G. 

3 

Its  worst  case  computational  complexity  is  0(n  ),  where  n  is  the  number 
of  nodes  in  G. 

In  [8],  K.  G.  Murty  and  C.  Perin  devleoped  an  efficient  lower 
bounding  strategy  for  the  minimum  objective  value  in  set  covering 
problems,  using  the  edge  covering  algorithm.  This  strategy  can  be  used 
to  develop  efficient  branch  and  bound  algorithms  for  the  set  covering 
problem. 

7.  THE  1-MATCHING/ COVERING  PROBLEM. 

Let  G  •  (N,  A)  be  an  undirected  network  with  n  nodes  and  m  edges, 
with  c  ■  (c^j )  as  the  edge  cost  vector.  Let  W"  ,  AT,  ,  N°  be  a  given 
partition  of  the  set  of  nodes  N,  A  1-matching/ covering  with  respect 
to  this  partition  is  a  subset  of  edges  F^*A  satisfying  the  property 
that  for  each  1  e  N;  there  is  at  most  one  edge  in  F  incident  at  i  if 
1  e  exactly  one  edge  if  i  e  W",  and  at  least  one  edge  if  1  e  W«. 

The  1-matching/covering  problem  is  to  find  a  minimum  cost  1  matching/ 
covering.  In  this  Ph.D.  dissertation  [23]  and  in  [16], 


V 


4 

C.  Perin  extended  the  blossom  algorithm  of  [6]  into  an 

efficient  algorithm  for  solving  the  1-matching/covering  problem,  whose 

3 

worst  case  computational  complexity  is  0(n  ).  He  has  also  developed 

«• 

mother  blossom  algorithm  of  the  out-of-kilter  type  in  [23J*  for  solving 

the  1-matching/ covering  problem,  which  also  has  the  worst  case  computa- 

3 

tional  complexity  of  0(n  ).  Combining  these  two  algorithms,  he  has 
developed  very  effficient  routines  for  sensitivity  analysis  in  the 
1-matchlng/ covering  problem,  when  changes  occur  in  the  edge  cost  vector, 
or  the  node  partition,  or  when  new  edges  or  nodes  are  introduced,  or  when 
edges  or  nodes  are  eliminated. 

C.  Perin  has  prepared  a  software  package  in  Fortran  for  these 
algorithms  for  the  1-matching/covering  problem  which  has  been  tested 
extensively,  and  found  to  perform  very  well  on  even  large  networks. 

[17]  is  a  user's  guide  for  using  this  package,  which  contains  a  complete 
listing  in  an  appendix. 

In  [15],  C.  Perin  describes  a  simple  proof  for  the  blossom 
constraints  characterizing  the  convex  hull  of  the  1-matching/covering 
incidence  vectors,  which  does  not  depend  on  the  algorithm. 

8.  SPECIFIED  CARDINALITY  1 -MATCHING/ COVERING  PROBLEM. 

This  is  the  1-matching/coverlng  problem  with  an  additional  constraint 

on  its  cardinality  r.  It  is  required  to  solve  this  problem  treating  r 

as  an  Integer  valued  parameter.  In  [18],  K.  G.  Murty  and  C.  Perin 

describes  a  blossom  algorithm  for  this  problem,  whose  worst  case 

2 

computational  complexity  is  0(n  m  ). 


9.  P0LYN0MIA1LY  BOUNDED  ELLIPSOID  ALGORITHMS  FOR  CONVEX  QUADRATIC  PROGRAMMING 


In  [7],  S.  J.  Chung  and  K.  G.  Murty  have  developed  polynomi ally  bounded 

) 

ellipsoid  algorithms  for  convex  quadratic  programs,  or  equivalently  LCPs 
associated  with  positive  semldeflnlte  matrices,  in  which  all  the  data  is 
Integer. 

Y.  Fathi  [13]  did  a  computational  study  of  several  algorithms  for 
solving  LCPs  associated  with  positive  definite  symmetric  matrices.  In 
these  computational  tests  he  found  out  that  the  critical  index  algorithm 
discussed  in  [4]  is  the  most  efficient  for  solving  this  class  of  problems. 
The  ellipsoid  method  of  [7]  performed  very  poorly  in  comparison  with  the 
others.  These  results  confirm  the  observatln  that,  while  the  ellipsoid 
method  is  a  tremendous  theoretical  breakthrough,  it  is  not  very  practical 
for  solving  linear  or  convex  quadratic  programming  problems. 

10.  COMPLEXITY  OF  ADJACENCY  CHECKING. 

In  [19,  24],  S.  J.  Chung  has  developed  efficient  linear  time  algorithms 
for  checking  whether  two  given  1-matching /cove ring  incidence  vectors  are 
adjacent  on  the  convex  hull  of  all  1-matching/covering  incidence  vectors. 

In  [20,  24],  S.  J.  Chung  has  studied  the  problem  of  checking  the 
adjacency  of  two  given  0-1  feasible  solutions  on  the  convex  hull  of  all 
■>*1  feasible  solutions  of  the  following  problems  (1)  the  0-1  equality 
constrained  Knapsack  problem,  (il)  the  0-1  inequality  constrained  Knapsack 
problem,  (Hi)  the  general  0-1  integer  programming  problem,  (iv)  the 
general  set  covering  problem.  In  each  case,  he  has  developed  necessary 
and  sufficient  conditions  for  two  given  0-1  feasible  solutions  to  be  adja¬ 
cent  on  the  convex  hull  of  all  0-1  feasible  solutions  of  that  problem. 
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But  unfortunately.  It  turns  out  that  checking  whether  these  conditions  hold. 

Is  an  DP-complete  problem  In  each  case.  This  conclusively  establishes  that 
the  problem  of  adjacency  checking  of  two  given  extreme  points'  on  the  Integer 
polytope  associated  with  each  these  problems  in  NP-complete. 

In  [20,  24],  S.  J.  Chung  has  also  established  that  the  problem  of  checking 
whether  two  given  extreme  points  of  an  unbounded  convex  polyhedron  are  adjacent 
on  the  convex  hull  of  all  extreme  points  of  this  polyhedron,  is  itself  an 
NP-complete  problem.  A  special  case  of  this,  to  check  whether  two  given 
simple  chain  incidence  vectors  are  adjacent  on  the  convex  hull  of  all  simple 
chain  incidence  vectors,  is  also  an  NP-complete  problem. 

11.  PERTURBATION  RESULTS  FOR  LCPs  ASSOCIATED  WITH  POSITIVE  SEMIDEFINITE 

MATRICES. 

Consider  an  LCP  (q,  M)  in  which  M  is  positive  semideflnite.  Consider 
the  perturbed  LCP  (q,  M  +  e  I)  where  e  >  0.  In  [21],  A.  Gana  has  derived 
the  limiting  behavior  of  the  solution  of  the  perturbed  problem  as  e  converges 
to  zero. 

12.  COMPLEXITY  OF  PARAMETRIC  INTEGER  AND  NETWORK  FLOW  PROBLEM. 

In  [12],  P.  Carstensen  has  studied  the  behavior  of  the  number  of  slope 
changes  in  the  optimum  objective  function  of  parametric  cost  single  commodity 
network  flow  problem  and  parametric  cost  0-1  integer  programming  problems. 
There,  she  has  constructed  a  parametric  cost  minimum  cost  flow  problem  with 
2n+2  nodes,  O(n^)  arcs,  with  0(2n)  slope  changes  in  the  optimum  objective 
value  function,  for  each  positive  integer  n.  Similarly  she  has  constructed 
parametric  cost  0-1  integer  programming  problems  (even  Knapsack  problems) 
in  which  the  number  of  slope  changes  for  the  optimum  objective  function  value 


grows  as  2  where  n  is  the  number  of  0-1  variables  in  the  problem. 

clearly  establishes  that  in  both  these  classes  of  parametric  cost 
problems,  the  number  of  slope  changes  in  the  optimum  objective  value 
curve  grows  exponentially  with  the  size  of  the  problem. 

13.  SOME  NP-OOMPLETE  PROBLEMS  IN  LINEAR  PROGRAMMING. 

In  [10 J,  R.  Chandrasekharan,  S.  N.  Kabadi  and  K.  G.  Murty  have 
established  that  many  of  the  unresolved  problems  in  linear  and  quadratic 
programming  (for  example,  the  problem  of  checking  whether  an  LP  with 
integer  data  is  degenerate,  whether  there  exists  a  basic  feasible  solution 
for  a  given  LP  with  a  specified  objective  value,  whether  a  given  square 
matrix  has  a  singular  principal  submatrix,  the  bilinear  programing  problem, 
etc.)  are  either  NP-complete  or  NP-hard. 

14.  CONVEXITY  OF  PIECEWISE  LINEAR  FUNCTIONS. 

In  [14],  K.  G.  Murty  has  established  the  necessary  and  sufficient 
conditions  for  a  given  piecewise  linear  function  on  a  convex  polyhedron, 
to  be  convex.  Be  has  developed  an  algorithm  for  checking  whether  the 
given  piecewise  linear  function  is  convex,  which  requires  the  solution 
of  a  small  number  of  linear 


programs. 
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